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Diffusion-induced Ramsey narrowing that appears when atoms can leave the interaction region 
and repeatedly return without lost of coherence is investigated using strong collisions approximation. 
The effective diffusion equation is obtained and solved for low- dimensional model configurations and 
three-dimensional real one. 
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I. INTRODUCTION 

The effect of diffusion- induced Ramsey narrowing has been considered in [1]. ft appears in the system of 
atoms moving through a laser beam in a buffer gas. The lifetime of the atomic coherence is limited and leads 
to the Lorentz line shape in assumption that atoms diffuse out of the laser beam and do not return. However, 
after spending some time in dark (outside the laser beam), atoms can return to the interaction region before 
losing the coherence (see fig. (1)). In some cases when decoherence effects are small diffusing atoms can spend 
a majority of their coherence time in the dark which leads to spectral narrowing of the center of the atomic 
resonance line shape in analogy to Ramsey spectroscopy [2] . 

In the present paper we propose the model for diffusion-induced Ramsey narrowing based on kinetic equation 
in approximation of strong collisions. The paper is organized as follows. In section II we discuss the basic kinetic 
equation in strong collision approximation. Further we obtain the effective diffusion equation in section III and 
solve it exactly in section IV for the simplest case unbounded configuration. Using initial approximation we 
obtain the similar diffusion equation for spatially bounded configuration in V. Short review of results is given 
in section VI. 
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Figure 1: Scheme of the experiment 
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II. BASIC EQUATIONS 



In the classical problem of migration of a particle in the gas the kinetic equation for the distribution func- 
tion p(t,v,t) has the form ([3], [4]): 

g + (tr.v)p=5, (i) 

where S is the collision integral that can be expressed in terms of p as follows: 

S = - J [A(v,if)p(f,v,t) - A(if,v)p(r,if,t)] dif . (2) 

Here A(v, if) means the numbers of transition v — > if during the unit of time due to collisions with the buffer 
gas. In this expression collisions between atoms are neglected. The initial condition is p(f, v, 0) = W(v)5(r), 
where W(v) is Maxwell distribution. 
The collision integral can be written as 

S = -vp{r, v,t) + J [- A(v, v)p(f, if, t)] dif , (3) 

where 

v{v) = j A(v, if) dv 1 (4) 

denote the effective collision frequency. 

There are two simple models that allow one to investigate the kinetic equation (1) with the collision integral 
(3). The first one based on the approximation of weak collisions (heavy atoms scatter on light particles) and 
leads to diffusion approximation with Fokker-Planck equation for the function p. In this model the velocity of 
the atom changes essentially after great number of collisions. 

The second model based on the approximation of strong collisions (light atoms scatter on heavy particles). 
In this model the function A(v' ,v) is independent of if, this means that the particle velocity after collision v is 
independent of its velocity before the collision. The form of A(v) can be determined from the condition S = 
in statical equilibrium. In this case v = const, and 

A(v) = vW(v) . 

This expression means that an arbitrary initial distribution turn to Maxwell after the lapse of time \ jv (or, 
cquivalently, the atom forgets its initial velocity as time 1/v passes). 

Finally, the kinetic equation for the problem of migration of the atom in gas in approximation of strong 
collisions has the form: 



— + (v ■ V)p = —vp + vW(v) j p(f,v,t)dv. 



The same approach can be used (see [4]) for non-diagonal elements of density matrix. For two- level system 
we have: 



dpu 
dt 



+ (v ■ V)pi2 = ~{v + 7 + iAcj)p 12 + vW{v) j pi 2 (r,v,t) dv + q(r,v), 



where 7 describes a spontaneous decay, q(f,v) is coherent excitation (with laser beam), Auj is two-photon 
frequency detuning. 

We will consider the stationary case. Denoting p{r,v) — pi2(r,v) we can write: 



(v ■ V)p(r, v) = —(y + 7 + iAuj)p(r, v) + W(v) 
where 

W{v) = W^e-^K , A(f) = \ Q e-^' a 
and Wq is the normalization constant for the Maxwell distribution. 



A(r) + v J dv p(r, v) 



(5) 
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Let us denote 



N(r) — J dv p{r, v) , a = v + 7 + iAu; , a = 7 + * Aw . 
The complex signal is denned as: 

S(Aw) = J J drdvX(r)p(r, v) = J N(f)X{f) dr. (6) 

We assume that 

a v 1 

Aw,7<^, f > — = — , 
a t 

where r a = ^- — is characteristic radiation zone time of flight. The last relations can be rewritten as r a t v = £ 
(collision time). 

For the atomic ensemble the atomic resonance line shape is determined by weighted average of the line shapes 
from different Ramsey sequences (with alternating interactions with the laser beam and motion in the dark). 



III. EFFECTIVE DIFFUSION EQUATION 

According to (II) the complex signal S(Aoj) can be expressed in terms of zeroth moment N(f) of the distribu- 
tion p(r) with respect to v. Let us obtain the equations for higher moments N^(r). For the sake of simplicity 
first we consider f-dimensional case, where 

00 

— CO 

Multiplying the equation 

v—p{x, v) = -ap(x, v) + W(v) [X(x) + ^V (0) (x)] 



by v k and integrating over v we get 



dx 



aiVW (x) + (v k ) [X(x) + vNM {x)] , (v k ) = J v W(v) dv . 



Since (v k ) = for odd values of k, we can write the following equations for odd and even orders: 

J7\r(2fc+1) J N (2k+2) 

= -aN^ + {v 2k } f X(x) + vN {0) {x) l t = - aN (^+i) 

dx dx 

or, equivalently equations for even orders only: 

1 d 2 7V(2fe+2) 



a 



. - aM 2k *> - (v 2k ) [X(x) + vN<M (x)] , (7) 



In these sequence of equations one can substitute the next equation to the previous and obtain the expression 
like the following: 

*t(0) rw n „foi/ m ld 2 N^ (v 2 ) d 2 rw , „ m . ,, 1 d 4 N^ 
aNi ) - ^ + ^ ^ = a— = c^"dx^^ + vNl0) W + ^ — - • ■ ' < 8 > 

Performing the substitution k times we obtain: 

(v 2 ) d 2 (v 2k ) d 2k \ r f0) , 1 d 2k+2 N^ 2k + 2 ^ 



« w<0, = ( 1 + ^^ + - + ^^) [AW+ ^° ,(l)1 + 



Q,2fe+1 ^ x 2k+2 
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For k — > oo this expression can be written symbolically: 



^ (0) =E^^[M*) + ^ (0) (*)] - F[A(x)+,iV(°)(x)] , (9) 

fe=0 



where F is a formal linear operator defined as 



oo 

J dvW(v 



v d x 



' a dx 

For higher dimensions the result is similar but more complicated because the moments became tensors: 

N it ti k (?) = J Vn-.. v lk p(r, v) dv 
In 3-dimensional case the same procedure leads to 

Hence, 

^V (0) = + -^(viVjjViVj + ^(viVjVkVijViVjVkVi + . . (A + vN<®) . 



Since 



(10) 



we can write finally 

#= ( i+ § A+ S Ai+ -)( A+ ^- 

The expressions (8) and (10) are asymptotic expansions and for sufficiently large v (and a) we can leave only 
first terms. The terms with second-order derivatives lead to the diffusion equation. It is easily to find that these 
terms have the same coefficients for all dimensions: 



v 2 

aN=(X + vN) + -^A(A^iV) + . . 
2a z 



or one can write up to the terms 1/v: 



vvl 



a N{f) = ^ AN(r) + A(f) (11) 
This equation can be treated as diffusion equation with absorption factor a = j + iAu; and diffusion coefficient 

D=^l~ V A (12) 
2a 2 2v V 1 

(for large v). 



Using the formula a — ^JDtd let us define diffusion characteristic time 

a 2 2a 2 

More suitable form of the diffusion equation (11) is 



TD = d " -zr v (13) 



AN(r)- /3 2 N(r) = -^A(r), /3 2 - . (14) 



a 

This equation can be solved using Green functions techniques. 
In the present system there are four characteristic times: 

a 1 1 T a 2 
T a = , ^7 = -, T v =-, T D ~ = VT a , 

v 7 ^ ^ 

According to initial assumption r a > t v (great number of collisions during the flight time) Also we assume 

T a < T 7 , T a <T D . 



4 



IV. EXACT SOLUTION FOR INFINITE REGION 



For spatially infinite region the equation (5) can be solved using Fourier transformation with respect to r: 

p{k,v)=[\{k) + vN{k)] — V , where p(k,v) = [ p(f,v) e^' p dr . 

a + ik ■ v J 

Since N(k) = J p(k,v)dv, we get after integration over v 

N(k)=[X(k) + uN(k)]F(k), where F(je) = f W ( v \ dv ; (15) 

J a + ik ■ v 

so that we can express the function N(k) as 

m = mP ® . as) 

l-vF(k) 

The function F(k) is the Voigt profile and can be presented as follows 

m = m, r . u {-fM _ a( ) d( (17) 

J a + ik-v J V 4 / 



This expression is valid for all spacial dimensions. Note that (17) can be written in terms of the error function 

a 
\k\vo 



namely F(k) = ^ Ae A crfc(A), where A = ju^- For sufficiently small k there exists the following asymptotic 



expansion (see [5]): 



1 



1 / kvo \ 2 3 / fc«o x 4 



a J 4 \ a 



where 



kv 




kv 


a 




V 



<1. (18) 



dkN(k)X{k), (19) 



Let us find the signal Soo . Taking into account the properties of the Fourier transformation we can write: 

S co (A^) = J dxN(r)X(r) = ^ J « 

where n denotes dimension of the space (1, 2 or 3). Substituting (16) we get 

X 2 (k)F(k) 



1 ' (27r)"i l-vF(k) 



;F(k) 

Due to the factor A 2 (fc) the value of S(Aw) is determined by small k, so that we can replace F(k) by its 
expansion (18): 

_ . . , 1 f X 2 (k)dk 

S TO Aw = — — — / V > k2v2 (20) 

2acto 

The final result depends on space dimension, we will consider partial cases n = 1 and n = 2 (the case n — 3 is 
unreachable on experiment if A(r) = A e~ r l a and is singular if A(r) = Aoe - ^ +v ^ a ). 



A. One dimension 

In this case 

W(v) = e- v2 / v « , X(x) = Aoe-" 2 /" 2 



•5 



bmce A(/c) = v ^a\ e- k ' a / 4 , we obtain from (20): 



+00 



Soo(Aw) 



X 2 a 2 f e~ k a ' 2 dk _ ira\ 2 Q A 2 



-00 2oa 



Ae A erfc(A), A 



2 aaoa 2 aaga 2 



2(v 2 ) ' 



(21) 



Since (v 2 ) = v 2 /2. Using the notations introduced above and (21) we can write Soo in the following form: 



S — 



7T aXlA 2 



2 a J x/T^K ' 


or, taking into account the fact that for v 3> 7, Aw the combination A 2 = aa^r 2 ~ aoru and 



S — 



aA r z 



^ [l + r/r D 



1/2 



(22) 



The right hand side can be interpreted as the superposition of profiles with the weight factor e Q ° T . 
If 7T/J > 1 we can expand 



, 1 r 
1 -2^ + '- 



^ \2 1 

aX — 

a 



1 1 
1-— ^— + ■ 



2 °«o + ^ 



This result is correct for small Aw (i. e. |Aw| <C 7) and jtd > 1, t < v. 
If 7Td ~ 1 we can expand directly in powers of Aw. For small Aw 



A 2 = a aT 2 ~ r D (7 + iAw) - t 2 Aw 2 , 



so that 



/To 



— e A crfc(A) = , /— e r ° 7 erfc(V^7) + — 



a 



7 

1 T D ( T D 



-77^ + ^(1 



1 



7 



1 /TD / T D 



— ^ \^--)+^xl^[^-^l-r 2 D -^- 2 



47 J 2 y 7 \ 7 



2r D 7 
3 

47 



e TD7 erfc(V^7) 
Aw 2 + . . . . 



iAw 



B. Two dimensions 



The Maxwell distribution function has the form 

W (v) = -L e-^K , X(x) = A e- ?2 / a2 

the Fourier transformation if A(r) is A(fc) = 7ra 2 A e~ fe2a2 / 4 , so that 



3oo(Aw) - — / 3 - — Ae Jiii (A J, ^ = — = , (23) 

2a 7 1 1 ^0 fc2 2a % {v^} 

2aao 



here Eii(. . . ) if the integral exponent of the first order [8]. 
The function can be rewritten as 



with the same meaning as (22). 
In the case 7775 3> 1 



-dor J x 

S^Aw) = / , (v 2 ) = v 2 . (24) 

TD 



r- /a s TrAga 2 1 
Soo (Aw) ~ 



2 «o + ^ 



G 



The combination vt 2 is the diffusion time for 1-dimcnsional problem. Note that in 2-dimensional configuration 
the real diffusion time defined as m — a 2 / D and expressed in terms of (v 2 ) is different (and can be obtained 
with replacing D — > 2D). 

If ~ 1 the expansion Soo (Aw) is 



— A 2 e A ' Ei(A 2 ) ~ T D e TD1 Ei(r D7 ) + — e r ° 7 [e TD7 + r D7 Ei(r z ,7)lzAw 
a 7 



Aw 2 



C. Interpretation of Diffusion Equation 

Using inverse Fourier transformation in (15) and the properties of the convolution we can write: 
N(r) = J[\{r') + vN{r')] F{r-r')dr' = j ' [\(r - f') + vN(r - f')] F(f')df', 

where 



F(r) 



(2tt)" 



F{k)e~ lk - r dk = W / exp -a£ 



+ 00 



^ V e 



Wo = 



i 



(V^o) n 



(25) 



(26) 



is transformation of F(k) (with the logarithmic singularity). Both of the convolutions in (25) can be approx- 
imated by an asymptotic expansion taking into account the peak if F(r) at r = 0), because its values are 
determined for f' close to r. Since 

N(f — r ') - N(f) + x' l V l N{r) + ^ x'^V.V 3 N{r) + ... 

we get 

J N(r-r')F(r')df" = N(r) j F(f") dr' + V t N(r) j x\F(r ') dr ' + 1 ViVjN(r) J x' l x' J F(f')dr" + ... 
Using the fact that F(r) depends only on r = |f| we can conclude that 

J N(r- f ') F(f") dr' = N(f) J F(f') dr' + 1 AN(f) J F(f')r' 2 dr' + ... 
Using the explicit form (26) one can obtain after simple calculations: 

J F{r') dr' = 1 , J F(f')r' 2 dr' = ^-J , (v 2 ) = j v 2 W{v) dv . 

Hence, 



N(r-r')F(r')dr' 



a 



(v 2 ) 

N(r) + ^AN(r) 



(27) 



The parameter | ^ | — ^ is assumed to be sufficiently small with respect to a. The same expansion is valid for 
convolution F(x) with X(x). Note that the convolution (25) is the same as (9). 

To prove that our approach is correct one can solve the diffusion equation and compare results with the 
previous calculations based on Fourier transformation. 



D. One dimension 



The Green function with the boundary conditions G(x, x')\ x =± 00 = is 

G(x,x') = ±e-^-*\ 
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so that 

+00 +00 



N(x) 



J_ I e -P\x-x '\f( x ')dx' = ^ [ e-^- x '\- x ' 2 / a2 dx' 
2/5 J 2a J 



4a 

+OO 



e"^ erfc - e) + crfc + e) 



where e = /3a/2. The integral SL> = / X(x)N(x) dx can be easily calculated and coincides with (21). 



E. Two dimensions 

The equation (14) can be solved using the Green function for R — ► 00 

G{r) = ±-K {(3r), 
Taking into account the boundary conditions we get 

N{r) = J dr'G(r-r')f{r') = ±- j dr' K a {(3\r - r'\)f{r') . 

The integral over angle variables can be calculated using the expansion 

00 

K ((3\f- r'\) = 7 (/3r<)^o(/3r>) + 2^ I n (/?r<)Jr n (/3r>) cosn(^ - ^) : 

71=1 

where r> and r< are the greater and smaller value from r and r' , 



\r — f'\ = \/r 2 + r' 2 — 2rr' cos{ip — ip') . 

Therefore 

N(r) = 2tt J dr' r7 (r < )X (r>) f(r') . (28) 

(2) 

The function Soo after simple transformations becomes 

S (2) = y rdrX ^ Ko ^ J d r'r'\(r')I (f3r'). (29) 


It can be shown (using the McDonald integral [6], 7.7.6 (37)) that this integral coincides with (23). 

V. DIFFUSION EQUATION IN FINITE REGION 

Other approach to diffusion equation is following: the equation (5) can be solved formally when the last 
term in right hand side is treated as inhomogeneous term (see [7]). The result is the same as (11). Namely, in 
1-dimensional case, where (5) is 



dp(x,v) . 
v p ^ ' ' + a oP (x, v) = W{v) 



\{x) + vN(x) 



, N(±R)=0. 



and —R^x^R The total solution reads 

p(x,v) = C(v)e- ax/v + 



X 

j EM [x( x ') + vN{x')] e - a ( x - x ')/ v dx' . 



Therefore the part of solution caused by inhomogencous term is 



oo +R 

»(,)=/** > ,)* = / i1 V) + V)] f (.-A 



F(x) was defined in (26). The homogeneous part of the solution that include C(v) vanishes according to initial 
assumptions due to the presence of the factors p(±i?, v) = — the coherence will be lost after collisions with 
walls x = ±R. 

The last integro-differential equation can be written in the form of diffusion equation. Taking into account 
the behavior of F(x) and expanding N(x) in series, one can obtain with additional assumption \aR/vo\ ^> 1 
(i. e. vR/v = vt r > 1): 



+R 

f 



N(x') F(x - x') dx' = 



1 



a 



N(x) 



a* 



N"(x) + ... 



this coincides with (27). Therefore, the diffusion equation is valid also for finite region when vtr 3> 1. 

Let us consider the solution of the diffusion equation in the form (14) for different dimensions. It appears 
that most effective for this purpose is the Green functions technique. 



A. One dimension 



The equation for N(x) is 



where 



N"(x) - (3 2 N{x) = -f{x) , N(-R) = N(R)=0. (30) 



f(x) = — A(r , (3 Z = - = . 

ao vv^ 2v{v z ) 

Taking into account the fact that both the functions N(x) and f(x) are even we can replace this problem by 

N"(x) -f3 2 N(x) = -f{x), N'iO) = NiR) = 0. (31) 

Let us find the Green G(x,x') function for this equation with the boundary conditions ^G(x, x')\ x=0 = and 
G(x, x')\ x= r = 0. It is easy to find that 



. cosh(/3x < ) smh(/?(i? - x>)) 

G(X, X ) = r , 

y ' (3cosh((3R) 

where x< and x < denote the greater and the smaller values from x and x' . 
The partial solution of inhomogeneous equation (31) is 



(32) 



Nix) = 



1 cosh/3a; 
/? cosh/3i? 



sinh/3(i? - x')fix') dx' 







sinh/3(x - x')fix') dx' 



Taking into account the boundary conditions we can conclude that this expression is a total solution of (31). 
Let us denote for simplicity 

X X 

C(x) = j e- y2 ^ 2 coshpydy, S(x) = J e^^ 2 sinh/3y dy . 



These integrals can be calculated in terms of error function: 

erf ^— — ej + erf ^— + s^j 



C{x)=^-a^ 
4 



C(oo) = ^ ae e 



Six)^a^ 
4 



erf (- - s) - erf (- + e) + 2 erf (e)l (oo) = ^ ae 5 ' erf (e) . 
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where e = 4p ~ A. Taking into account C'(x) = e x l a cosh/3a; and S'(x) = e x l a sinh/fc, the complex 
function can be rewritten as 



R a °J cosh/3 R y ' 



It is easy to show when R — > oo S^P — ► according to (21). Hence, current interpretation of the initial 
integral equation in terms of diffusion (31) equation is correct in our approximation. 

For i?> a one can compare with S^J . After simple calculations we can write 

S W = S «-^Mc 2 ( J R)[l-tanh(/3 J R)] 

For large x 

C(x) ~ C(oo) - ^ ae e2 erfc - e) 

therefore 

s£ } ~ S W + 27Fq2/3 e 2e2 e~ 2/3fl . 



B. Two dimensions 

The diffusion equation in the region r < R 

AN(f) - p 2 N(r) = -f(f) , N(r)\ r=R = 
can be solved using the Green function 

AG(r) - p 2 G(r) = -S(r- f ') , G(r)\ r=R = . 

where 

h(M T lo ,. v faos ^ I n (J3r) 



Ko ^~ J^jm) ^^KoiPR) ~ 2 ]T JJj^ InWr')K n (PR) cosn(^ - tf) 



All the terms with n^O vanish in the integrals over angles and the term with Ko((3\r — f'\) can be written as 
in the case with R — > oo. 

The expression with N(r) has the form 

,^^ / o(^<)[ / o(^)^o(^>) - IoiPr^KojpR)] ^ 



N(r) = f r'f(r') 
Jo 

r R R 

= K (Pr) J dr' r'I (pr')f(r') + I (/?r) J dr' r' K ((3r')f(r') - J (/?r) J dr' r> 7 (/3r')/(r') . (33) 



The function is 



S {2) = 2w J N(r)X(r)rdr. 



o 

Using the symmetry of the integrands Sr can be presented in more compact form: 

S ( 2 ) = 2/3^ K (fir)I (fiR) dr 

a J K (pR) 



o 
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where 



/2 / r>2 

xl (f3x)e- x2 / a2 dx , J(oo) = y exp ( ?-2L 



For large r (substituting r instead of x in the integrand): 

oo 

J(r) = J(oo) - y a;/o(/3a;)e" ;l:2/a2 ~ J(oo) 



Like to 1-dimensional case one can compare with Soo ■ The result is: 



;( 2 ) 



I Q (J3r) K a {(3r) 
I a (pR) K (fiR) 



Since K (x) ~ x , J (a:) 



/27TX 
OO 



for a; 3> 1 we can rewrite the last expression as 



Sg } = SP)+47rAg^- / V^Rsinh[/3(r-i?)le-'- 2 / a2 X(/3r)dr~7r 2 A2^ e - 2 ^X(oo) 

~ S P) - Aq— f2^e^- fi2 / a2 + tt 2 e" 2 ^) X(oo) 2 . 
an V / 



C. Tree dimensions 

Consider the case when N(r) vanishes on the cylinder of radius r = R bounded by the planes z = ±1. The 
diffusion equation remains the same: 

AN(r) - [3 2 N(r) = -f(r) , N(r)\ r=R = , N(r)\ z=±l = 

The Green axially symmetric function satisfies the equation 

AG(f, f') - (3 2 G(f, r ') = - 6{r ~ r '^ Z ~ Z ' ] , G(r,r% =R = 0, G(r, f ') | z=±J = . 

First consider the case R — > oo. Expanding G(r, r') in series of <p n {z) — sin ^ (I — z) (with <f> n (±l) = 0) we 
get 



^ OO 

= 7 ^ sin(fc„z) sin(/c„z')/ (A„r < )i : s:o(A„r>) , where k n 

ZTTl * — ; 



n=l 



21 



Pi = P Z + K. (35) 



Finally, the axially symmetrical solution of the diffusion equation is 

oo 00 

N{?)= J G(f,r)f(f') df' = ^Y1 S1 2n + T I rfrV M/Wr<)tfo(/Wr>)/(r') • 



n=l 



Note that the integral over r is the same that N(r) for 2-dimensional problem with (3 — > /32n+i (see (28)). 
The complex signal is 

oo 

S^=2TrJJdrdzrN(r)\(r) = ^E ( 2n +i)2 / drrA(r)JV(r) , 



or, denoting by S^(/3) the profile (29) for corresponding 2-dimensional case we get 

00 ^ J -^2^ (2n + i)2- 
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S(A«) 



— n=0 

— -n=1 

— n=2 





0,6 



Figure 2: Terms in power series comparison of with 



;( 2 ) c(3) 



The same procedure can be used for R < oo, the Green function is 



hence, 



rl = ^f; sin(fc„z) sin(fc„z') [7 (M #o(/3„r>) - I (/? n r>)Jr (/3nfl)] , (36) 

27T/ ^ lo(Pn-ti) 

oo fl 

S^ 3) = 2* J J drdzrN{r)\(r) = ^ £ (2n + 1)2 / dr rA(r)JV*(r) , 



again, denoting by S^(/3) the profile (34) for corresponding 2-dimensional case we obtain 

S (3) [5 1- 2 ^ V S R^n] 



n=0 



In this series only the first terms are significant for result (see fig. 2). 



VI. DISCUSSION 

The graphics of S(Aw) with fixed parameters 

v = 1000 , 7 = 1, v = 100 , a = 1 . 
are presented at the figures. Characteristic times are 

T a = 0.01 , T 7 = 1 , T D = 0.1 . 

The values of i? (numerically coincide with i?/a) arc R = 2.0 , 3.0 , 5.0. All the functions are normalized. One 
can see that in the signal in 2-dimensional case is broader than in 1-dimcnsional. 
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Figure 3: The function S(Au>) for R = oo and R < oo (dash), half- width as the function of R — 1-dimensional case 
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Figure 4: The function S(Au>) for R = oo and R < oo (dash), half- width as the function of R — 2- dimensional case 




Figure 5: The function S(Au) for R < oo with fixed R and I — 3-dimensional case 
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